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Ch ' Abstract 
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We report anisotropy of the upper critical field (B C 2) of an intercalated graphite 
superconductor YbC6 {T c = 6.5 K) determined from angular dependent magnetore- 
sistance measurements. Though the perpendicular coherence length is much longer 
than interlayer spacing, measured angular dependences of B C 2 are well fitted by the 
Lawrence-Doniach model or the Tinkham model, which are known to be applica- 
ble to quasi two-dimensional materials or thin films, rather than the effective mass 
model. This observation is similar to the measurements for the other intercalated 
graphite superconductor, CaCg, by Jobiliong et al. [E. Jobiliong, H.D. Zhou, J. A. 
Janik, Y.-J. Jo, L. Balicas, J.S. Brooks, C.R. Wiebe, Phys. Rev. B 76 (2007) 31 
052511]. A possible explanation for the unexpected applicability of these models is 
that our YbCg samples are synthesized as thin flakes in the host graphite. 
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Graphite intercalation compounds (GICs) have been intensively studied for 

many years. This is because the electronic and magnetic properties can be con- 
trolled to a large extent by changing intercalant species and their arrangements 
such as the stage number and the intralayer structure [TJ. Superconductivity 
is one of the most interesting properties among them. The superconducting 
transition temperature (T c ) of GICs had been long limited below 3 K under 
ambient pressure [1] . The intercalants here are usually alkali metals. Recently, 
it was discovered that CaC 6 and YbC 6 have much higher T c of 11.5 and 6.5 
K, respectively [2]. Measured properties of these new high-T c GICs are well 
described with the framework of BCS theory [3,4.5.6.7]. Although the phonon 
mechanism is supposed to be relevant for their Cooper pairing [8], detailed 
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microscopic mechanisms for the attractive interactions have not been fully 
understood yet. 

Since GICs are structurally anisotropic, the superconducting properties also 
show anisotropic features. For instance, the upper critical field (B c2 ) of KHgCs 
shows angular dependences which are well fitted by the effective mass (EM) 
model j9] . On the other hand, the anisotropy of B c2 in CaC6 is accounted for by 
the Lawrence-Doniach (LD) model or the Tinkham model [10], though these 
models are applicable mainly to quasi two-dimensional (2D) materials or thin 
films. Here we report measurements of angular dependent magnetoresistance 
in YbC6- The anisotropy of B c2 deduced from our data is well described by 
the LD model or the Tinkham model. We infer that this result comes from a 
flake-shaped distribution of YbC@ in the host graphite in our sample. 

The YbCe samples were synthesized from the host material, a highly oriented 
pyrolytic graphite (HOPG) [11], by the vapor transport method. A pellet of 
Yb and a piece of HOPG were vacuum-sealed in a glass ample, and heat- 
treated at 400 °C for 2 days. After cooling down, the glass ample was opened 
inside a glove box with N 2 flow, because YbCg is extremely reactive to air. In 
the case of samples for transport measurements, Au lead wires were attached 
on HOPG with Ag paste before the heating. The typical sample size is about 
5x2x0.3 mm 3 . 

The synthesized samples were first analyzed by X-ray diffraction (XRD) with 
CuKa radiation at room temperature. The sample surfaces for the XRD mea- 
surements were protected with adhesive tapes against reaction in air. Fig. 1 is 
a measured XRD spectrum for such a sample showing both YbC6 (001) peaks 
and host graphite ones. The interlayer spacing (d') deduced from this spec- 
trum is 0.456 nm which is almost the same as the value reported by Weller et 
al. [2]. 

DC magnetization measurements in magnetic fields (B) parallel to the ab- 
plane of YbCe were carried out with a SQUID magnetometer [12]. The inset 
of Fig. 2 shows temperature dependences of magnetization (M) measured at 
B = 1 niT by field cooling (FC, closed circles) and zero field cooling (ZFC, 
open circles). It shows T c = 6.5 K, and an estimated superconducting volume 
fraction is a few percent. Field dependences of M measured at various fixed 
temperatures are shown in Fig. 2. In this plot, linear background signals, which 
are at most 6xl0 -3 emu at B — 1 T, have already been subtracted from the 
raw data. The origin of this background signal is not known at present. At 
each temperature, the upper critical field in the parallel direction (-B C 2||) is 
determined as a field where magnetization returns to zero as indicated by the 
arrows in the figure. 

Resistance measurements across the afe-plane were carried out using the stan- 



dard four-probe method [13] down to T = 2 K. A measurement current was 
fixed at 1 mA. Sample surfaces for the transport measurements were protected 
with epoxy resin against reaction in air. Fig. 3 shows the temperature depen- 
dence of resistance (R) in zero magnetic field. The T c value (= 6.47 K) here 
is defined as R(T C ) = \ R(T = 6.6 K) (see the inset of Fig. 3). This is con- 
sistent with the T c value determined from the magnetization measurements. 
The transition width corresponding to a resistance change from 0.1 to 0.9i?(T 
= 6.6 K) is 0.1 K. Note that R decreases to zero at T < 6.4 K, implying 
that superconducting paths connect the two voltage leads in spite of the small 
volume fraction of YbCg in the host. 

Field dependences of resistance measured at various fixed temperatures are 
shown in Fig. 4. The field is applied parallel to the afr-plane. At each tem- 
perature, the resistance increases rather sharply above a certain field which 
corresponds to B c2 \\ as will be discussed later. With further increasing field, it 
approaches gradually to a universal curve (R n (B)) in the normal state whose 
apparent field dependence is governed by the host material. Note that the 
magnetoresistance of graphite is very large 



B C 2\\ values determined in this way are plotted as a function of temperature 
by the closed circles in Fig. 5. B C 2\\ varies linearly with T in the measured 
temperature range (2 K < T < T c ). There are in good agreement with those 
determined from the magnetization measurements shown by the crosses in 
the figure. This consistency suggests appropriateness of the two independent 
determinations of B c2 \\. The open circles are fields where R{B) = | R n (B). 
This is not an appropriate definition for B c2 \\, since they are substantially 
(by ~ 65 %) larger than the B c2 \\ values determined from the magnetization 
measurements. 

Next, we show results of resistance measurements in magnetic fields of various 
directions in Fig. 6. We define 9 as an angle between the field direction and 
the c-axis of YbC 6 . As 9 decreases from 90° to 0°, the field dependence of 
R n (B) becomes steeper reflecting the anisotropic magnetoresistance of host 
graphite [14]. From these data, temperature dependences of B c2 in various 
field directions are obtained (Fig. 7). The T-linear dependence of B c2 holds 
for any 9. Anisotropy parameter (e = B c2 \\/B c2 ±) is calculated to be 2.1 for 
T = 2 K. Angular dependences of B c2 at various fixed temperatures are also 
obtained from the data in Fig. 6 and plotted in Fig. 8. At any temperature, 
B c2 is anisotropic and a cusp is observed around 9 = 90°. 

We now compare the measured properties of YbCg with three different the- 
oretical models for anisotropic superconductivity in layered materials or thin 
films. The first one is the EM model [15], which introduces the anisotropy of 
the effective mass into the GL equations. This model deals superconductors as 
anisotropic continua, and is relevant when the coherence length perpendicular 



to layers is much longer than the interlayer spacing. In this model, the angular 
dependence of B c2 is given by 



' B c2 (9)sm6 \ 2 ( B c2 {0) cos O V 

B C 2\\ J \ B c2 ± J 



1. 



where B c2 ± is the upper critical field in perpendicular direction. Also, relations 
between B c2 and the coherence length (£) are given by 



and 



B -< T >=s|b = diw( 1 -I) iTBsT ' ] (2) 

Bc2|(T) = 2^||(T)a(T) = 2^,(0)^(0) I 1 " 5;) (T * Tc) ' (3) 

Here $o — ^/2e is the flux quantum, and £j_ and £11 are coherence lengths 
perpendicular and parallel to layers, respectively. From Eq. (J2J), Eq. ([3} and 
the data in Fig. 7, £j_ and £y are estimated as 

£i(0) = 17 nm, (4) 

£|,(0) =37nm. (5) 

This £j_(0) value is much larger than the interlayer spacing d! = 0.456 nm, 
which does not exclude the applicability of the EM model. However, as shown 
in Fig. 8, fitting of the angular dependences of B c2 to Eq. (jTJ is poor (dashed 
lines). Particularly, this model does not represent the cusps around 9 = 90°. 

There are two other models which may be applicable to our data. They are 
the LD model [16] and the Tinkham model [T7] . The LD model describes 
anisotropic superconductors as stacked 2D superconducting layers connected 
by Josephson coupling. On the other hand, the Tinkham model is to describe 
superconducting thin films. In both models, the angular dependence of B c2 is 
given by 



B c2 {6) sin6 



+ [ Z7 ) =L (6) 



b c2 (6) cose y = 1 

B c2 \\ 

These two models provide much better fittings for the measured angular de- 
pendences in Fig. 8 (solid lines). This suggests the thin film nature of our 
sample. Since £j_ ^> d', it would be hard to apply the LD model to the present 
results. Instead, YbC^ might be synthesized as thin flakes in the host acting as 
independent superconducting thin films. This is conceivable for the following 
reason. The host HOPG contains high-density stacking faults, and its effective 
thickness is estimated to be around 40 layers (~ 13 nm) [18]. So, it is likely 
that the thickness of synthesized YbCg is also restricted by this layer number. 

Although the Tinkham model is originally developed for the case of d < £j_ 
[T7] , it is known to be a fairly good approximation even for the case of d ~ £j_ 



where e ~ 2 unless d ^> £j_ [19] . Therefore, it seems reasonable that the 
angular dependences of B c2 obey Eq. (jSJ), while e (= 2.1) is not so large. The 
same angular dependences of B c2 were also reported by Jobiliong et al. [TU] for 
another high-T c GIC superconductor CaCg (T c = 11.5 K), where e = 3.7 ~ 4.6. 

With the Tinkham model, other measured properties of our sample can be 
explained. First, a possible thickness variation in the superconducting paths 
would broaden the resistive transition in magnetic fields, since B c2 increases 
with decreasing thickness in thin film superconductors [20] . Next, the fact that 
the B c2 values measured in this work are larger than those by Weller et al. [2] 
for 9 = 90° (see Fig. 7) may be also due to the thinner effective thickness of 
our sample. However, such a difference is not seen at T > 4 K for 9 = 0°, for 
which we don't have reasonable explanation at this moment. 

There remain a few points which are not accounted for by the Tinkham model. 
For instance, the measured B c2 at 9 = 90° shows T-linear dependence rather 
than y/T c — T dependence 

B c2 (T) = v/12^ oc ^T c -T (7) 

predicted by this model. Here, d is film thickness. The same disagreement has 
been reported by Jobiliong et al. [10] even for a bulk GIC superconductor, 
i.e., CaCe synthesized from HOPG by the direct reaction method [10] . On the 
other hand, for YbC6, the temperature dependence of B c2 at 9 = 90° measured 
by Weller et al. [2] is not linear with T nor y/T c — T. In order to solve this 
puzzle, future synthesis of YbC 6 from single crystal graphite in well-defined 
form such as bulk is highly desirable. 

In conclusion, we synthesized YbC 6 , a superconducting GIC with T c = 6.5 K, 
by the vapor transport method. The measured angular dependences of B c2 at 
various field directions are well fitted by the Lawrence-Doniach (LD) model 
or the Tinkham model. This indicates the thin film nature of sample, which 
might originate from the thin film nature of host HOPG with high-density 
stacking faults. This scenario including the currently unexplained temperature 
dependence of B c2 will be tested by synthesizing sample with better defined 
form in future. 
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Fig. 1. The XRD spectrum of a synthesized YbCg sample measured with CuKa 
radiation at room temperature. Only (001) peaks are observed because of the c-axis 
orientation of both host HOPG and YbCg- 
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Fig. 2. The field dependence of magnetization in YbCg- The field direction is parallel 
to the afc-plane. Linear background signals have already been subtracted from the 
raw data. -B C 2|| values are indicated by the arrows. Inset: the temperature depen- 
dence of magnetization measured at B = 1 mT (B \\ arj-plane) with field cooling 
(FC) and zero field cooling (ZFC). 
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Fig. 3. The temperature dependence of resistance for YbCg in zero magnetic field. 
Inset: the detailed resistance change near T c (= 6.47 K). 
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Fig. 4. The field dependence of resistance in YbCg- The field direction is parallel to 
the afr-plane. The arrows indicate B c2 \\ values where resistance becomes finite with 
increasing field. 
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Fig. 5. The B c2 \\ values of YbCg determined from the transport measurements (cir- 
cles) and the magnetization measurements (crosses). The open circles are fields 
defined as midpoints of the resistive transition, which are substantially larger than 
the B c2 \\ values (see text). 
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Fig. 6. The magnetic field dependences of resistance for YbCg in various field direc- 
tions. 6 is denned as an angle between the magnetic field direction and the c-axis 
of YbC 6 . (a) 9 = 90°; (b) 9 = 60°; (c) 9 = 30°; (d) 9 = 0°. 
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Fig. 7. The temperature dependences of B& for YbC6 in various field directions. 
The T- linear dependence of B C 2 holds at any 6. The dashed lines are from Weller 
et al. [2]. 



14 




90 120 

[cleg] 



180 



Fig. 8. The angular dependence of B C 2 in YbCg at various fixed temperatures. 
Fittings with the EM model and the LD or Tinkham model are shown by the 
dashed lines and solid lines, respectively. 
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